Cardioviral RNA structure logo analysis: entropy, correlations, and prediction

In recent years, there has been an increased number of sequenced RNAs leading to the development of new RNA databases. Thus, predicting RNA structure from multiple alignments is an important issue to understand its function. Since RNA secondary structures are often conserved in evolution, developing methods to identify covariate sites in an alignment can be essential for discovering structural elements. Structure Logo is a technique established on the basis of entropy and mutual information measured to analyze RNA sequences from an alignment. We proposed an efficient Structure Logo approach to analyze conservations and correlations in a set of Cardioviral RNA sequences. The entropy and mutual information content were measured to examine the conservations and correlations, respectively. The conserved secondary structure motifs were predicted on the basis of the conservation and correlation analyses. Our predictive motifs were similar to the ones observed in the viral RNA structure database, and the correlations between bases also corresponded to the secondary structure in the database.

The topological configuration and conformational analysis of mRNA in translation

The theoretical model construction of mRNA hairpin structure and single-stranded structure as well as the simulation studies on RNA structure determined by the X-ray crystal diffraction and nuclear magnetic resonance revealed that in translation, after mRNA being unfolded into single-stranded structure, its topological configuration was closely correlative with the original hairpin structure. The conformational features of single-stranded mRNA appeared as helical regions alternating with curly regions to different extents, which might exert the influence on the folding of nascent polypeptide by various regulating effects including different translational rates.

Statistical correlation between protein secondary structure and messenger RNA stem-loop structure

Anew integrated sequence-structure database, called IADE (Integrated ASTRAL-DSSP-EMBL), incorporating matching mRNA sequence, amino acid sequence, and protein secondary structural data, is constructed. It includes 648 protein domains. Based on the IADE database, we studied the relation between RNA stem-loop frequencies and protein secondary structure. It was found that the alpha-helices and beta-strands on proteins tend to be preferably "coded" by mRNA stem region, while the coils on proteins tend to be preferably "coded" by mRNA loop region. These tendencies are more obvious if we observe the structural words (SWs). An SW is defined by a four-amino-acid-fragment that shows the pronounced secondary structural (alpha-helix or beta-strand) propensity. It is demonstrated that the deduced correlation between protein and mRNA structure can hardly be explained as the stochastic fluctuation effect. (C) 2003 Wiley Periodicals, Inc.

Self-assembly of single-stranded RNA on carbon nanotube: Polyadenylic acid to form a duplex structure

All messenger-RNA (mRNA) molecules in eukaryotic cells have a polyadenylic acid [poly (rA)] tail at the 3'-end and human poly (rA) polymerase (PAP) has been considered as a tumor-specific target. A ligand that is capable of recognizing and binding to the poly(M) tail of mRNA might interfere with the full processing of mRNA by PAP and can be a potential therapeutic agent. We report here for the first time that single-walled carbon nanotubes (SWNTs) can cause single-stranded poly (M) to self-structure and form a duplex structure, which is studied by UV melting, atomic force microscopy, circular dichroism spectroscopy, and NMR spectrometry.

Stable RNA hairpins in 88 coding regions of human mRNA

RNA hairpins containing UNCG, GNRA, CUUG (N = A, U, C or G, R = G or A) loops are unusually thermodynamic stable and conserved structures. The structural features of these hairpin loops are very special, and they play very important roles in vivo. They are prevalent in rRNA, catalytic RNA and non-coding mRNA. However, the 5' C(UUCG)G 3' hairpin is not found in the folding structure of 88 human mRNA coding regions. It is also different from rRNA in that there is no preference for certain sequences among tetraloops in these 88 mRNA folding structures.

Transcriptome of embryonic and neonatal mouse cortex by high-throughput RNA sequencing

Brain structure and function experience dramatic changes from embryonic to postnatal development. Microarray analyses have detected differential gene expression at different stages and in disease models, but gene expression information during early brain development is limited. We have generated >27 million reads to identify mRNAs from the mouse cortex for>16,000 genes at either embryonic day 18 (E18) or postnatal day 7 (P7), a period of significant synapto-genesis for neural circuit formation. In addition, we devised strategies to detect alternative splice forms and uncovered more splice variants. We observed differential expression of 3,758 genes between the 2 stages, many with known functions or predicted to be important for neural development. Neurogenesis-related genes, such as those encoding Sox4, Sox11, and zinc-finger proteins, were more highly expressed at E18 than at P7. In contrast, the genes encoding synaptic proteins such as synaptotagmin, complexin 2, and syntaxin were up-regulated from E18 to P7. We also found that several neurological disorder-related genes were highly expressed at E18. Our transcriptome analysis may serve as a blueprint for gene expression pattern and provide functional clues of previously unknown genes and disease-related genes during early brain development.

Functional domains and the antiviral effect of the double-stranded RNA-dependent protein kinase PKR from Paralichthys olivaceus

The double-stranded RNA (dsRNA)-dependent protein kinase PKR is thought to mediate a conserved antiviral pathway by inhibiting viral protein synthesis via the phosphorylation of the alpha subunit of eukaryotic initiation factor 2 (eIF2 alpha). However, little is known about the data related to the lower vertebrates, including fish. Recently, the identification of PKR-like, or PKZ, has addressed the question of whether there is an orthologous PKR in fish. Here, we identify the first fish PKR gene from the Japanese flounder Paralichthys olivaceus (PoPKR). PoPKR encodes a protein that shows a conserved structure that is characteristic of mammalian PKRs, having both the N-terminal region for dsRNA binding and the C-terminal region for the inhibition of protein translation. The catalytic activity of PoPKR is further evidence that it is required for protein translation inhibition in vitro. PoPKR is constitutively transcribed at low levels and is highly induced after virus infection. Strikingly, PoPKR overexpression increases eIF2 alpha phosphorylation and inhibits the replication of Scophthalmus maximus rhabdovirus (SMRV) in flounder embryonic cells, whereas phosphorylation and antiviral effects are impaired in transfected cells expressing the catalytically inactive PKR-K421R variant, indicating that PoPKR inhibits virus replication by phosphorylating substrate eIF2 alpha. The interaction between PoPKR and eIF2 alpha is demonstrated by coimmunoprecipitation assays, and the transfection of PoPKR-specific short interfering RNA further reveals that the enhanced eIF2 alpha phosphorylation is catalyzed by PoPKR during SMRV infection. The current data provide significant evidence for the existence of a PKR-mediated antiviral pathway in fish and reveal considerable conservation in the functional domains and the antiviral effect of PKR proteins between fish and mammals.

Variation patterns of the mitochondrial 16S rRNA gene with secondary structure constraints and their application to phylogeny of cyprinine fishes (Teleostei : Cypriniformes)

The mitochondrial 16S ribosomal RNA (rRNA) gene sequences from 93 cyprinid fishes were examined to reconstruct the phylogenetic relationships within the diverse and economically important subfamily Cyprininae. Within the subfamily a biased nucleotide composition (A > T, C > G) was observed in the loop regions of the gene, and in stem regions apparent selective pressures of base pairing showed a bias in favor of G over C and T over A. The bias may be associated with transition-transversion bias. Rates of nucleotide substitution were lower in stems than in loops. Analysis of compensatory substitutions across these taxa demonstrates 68% covariation in the gene and a logical weighting factor to account for dependence in mutations for phylogenetic inference should be 0.66. Comparisons of varied stem-loop weighting schemes indicate that the down-weightings for stem regions could improve the phylogenetic analysis and the degree of non-independence of stem substitutions was not as important as expected. Bayesian inference under four models of nucleotide substitution indicated that likelihood-based phylogenetic analyses were more effective in improving the phylogenetic performance than was weighted parsimony analysis. In Bayesian analyses, the resolution of phylogenies under the 16-state models for paired regions, incorporating GTR + G + I models for unpaired regions was better than those under other models. The subfamily Cyprininae was resolved as a monophyletic group, as well as tribe Labein and several genera. However, the monophyly of the currently recognized tribes, such as Schizothoracin, Barbin, Cyprinion + Onychostoma lineages, and some genera was rejected. Furthermore, comparisons of the parsimony and Bayesian analyses and results of variable length bootstrap analysis indicates that the mitochondrial 16S rRNA gene should contain important character variation to recover well-supported phylogeny of cyprinid taxa whose divergences occurred within the recent 8 MY, but could not provide resolution power for deep phylogenies spanning 10-19 MYA. (c) 2008 Published by Elsevier Inc.

Mutations stabilize small subunit ribosomal RNA in desiccation-tolerant cyanobacteria Nostoc

The ribosomal RNA molecule is an ideal model for evaluating the stability of a gene product under desiccation stress. We isolated 8 Nostoc strains that had the capacity to withstand desiccation in habitats and sequenced their 16S rRNA genes. The stabilities of 16S rRNAs secondary structures, indicated by free energy change of folding, were compared among Nostoc and other related species. The results suggested that 163 rRNA secondary structures of the desiccation-tolerant Nostoc strains were more stable than that of planktonic Nostocaceae species. The stabilizing mutations were divided into two categories: (1) those causing GC to replace other types of base pairs in stems and (2) those causing extension of stems. By mapping stabilizing mutations onto the Nostoc phylogenetic tree based on 16S rRNA gene, it was shown that most of stabilizing mutations had evolved during adaptive radiation among Nostoc spp. The evolution of 16S rRNA along the Nostoc lineage is suggested to be selectively advantageous under desiccation stress.

The molecular characterization of RNA segment S9 of grass carp hemorrhage virus (GCHV), an aquareovirus

Although reovirus infection is one of the major virus diseases of grass carp in China, the available knowledge on the structure and function of genes and proteins of the virus is limited. The complete sequence of the S9 genome segment of grass carp hemorrhage virus (GCHV) was determined. The segment consists of 1130 nucleotides and has a large open reading frame (ORF) encoding a protein of 352 amino acids with predicted molecular mass of 37.7 kDa. Amino acid sequence comparison revealed that the deduced protein encoded by GCHV S9 is closely related to the sigma NS proteins of mammalian reovirus (MRV) and avian reovirus (ARV). Secondary structure analysis displayed that the form of alpha -helices (40.1%) and beta -sheets (49.4%) are the richest two contents in the protein encoded by S9, and this protein is predicted to be a nonstructural protein. (C) 2001 Elsevier Science B.V. All rights reserved.

Evaluation of different culture conditions for high-level soluble expression of human cyclin A(2) with pET vector in BL21 (DE3) and spectroscopic characterization of its inclusion body structure

In this paper, we evaluated various parameters of culture condition affecting high-level soluble expression of human cyclin A, in Escherichia coli BL21(DE3), and demonstrated that the highest protein yield was obtained using TB(no glycerol) + 0.5% glucose medium at 25 degrees C. By single immobilized metal ion affinity chromatography, we got highly purified human cyclin A(2) with a yield ranged from 20 to 30 mg/L. By amyloid-diagnostic dye ThT binding and Fourier transform infrared spectroscopy, we observed a significant decrease in alpha-helix content and an increase in beta-sheet structure in cyclin A(2) inclusion body in comparison to its native protein, and confirmed the resemblance of the internal organization of cyclin A(2) inclusion body and amyloid fibrils.

Cloning and structure analysis of hydrogenase gene from Chlamydomonas reinhardtii SE

The cDNA of Chlamydomonas reinhardtii SE encoding hydrogenase (HydA2) was obtained from the total RNA of C reinhardtii SE by RT-PCR. The DNA of hydrogenase was amplified by PCR from the genomic DNA of C reinhardtii SE. The cDNA and DNA of hydrogenase were sequenced, respectively. The structure of hydrogenase gene was analyzed by biology software. The open reading frame predicts that the hydrogenase is composed of 3584 bp encoding 505 amino acids in length with a predicted M.W. of 53.69 kDa. Ten exons (including 1518 bp) and nine introns (including 2066 bp) have been found in the hydrogenase, and there were two potential N-glycosylate sites, eight protein kinase C phosphorylation site, eight casein kinase H phosphorylation site and one sulphorylation in the sequence. The theory pI was 6.15. Total number of negatively charged residues (Asp + Glu) and positively charged residues (Arg + Lys) were 55 and 61, respectively. (c) 2005 Elsevier Ltd. All rights reserved.